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> ■ ABSTRACT 

V^Q ' The power spectrum of weak lensing shear caused by large-scale structure is an emerg- 

CNJ . ing tool for precision cosmology, in particular for measuring the effects of dark energy 

on the growth of structure at low redshift. One potential source of systematic error is 
O'N intrinsic alignments of ellipticities of neighbouring galaxies (II correlation) that could 

mimic the correlations due to lensing. A related possibility pointed out by Hirata 
and Seljak (2004) is correlation between the intrinsic ellipticities of galaxies and the 
density held responsible for gravitational lensing shear (GI correlation). We present 
• constraints on both the II and GI correlations using 265 908 spectroscopic galaxies 

Q-f from the Sloan Digital Sky Survey (SDSS), and using galaxies as tracers of the mass in 

the case of the GI analysis. The availability of redshifts in the SDSS allows us to select 
£5 ' galaxies at small radial separations, which both reduces noise in the intrinsic align- 

^ | ment measurement and suppresses galaxy-galaxy lensing (which otherwise swamps 

C3 • the GI correlation). While we find no detection of the II correlation, our results are 

nonetheless statistically consistent with recent detections found using the SuperCOS- 
MOS survey. Extrapolation of these limits to cosmic shear surveys at z ~ 1 suggests 
■ that the II correlation is unlikely to have been a significant source of error for current 

' measurements of ag with ~ 10 per cent accuracy, but may still be an issue for future 

surveys with projected statistical errors below the one per cent level unless eliminated 
using photometric redshifts. In contrast, we have a clear detection of GI correlation in 
galaxies brighter than L* that persists to the largest scales probed (60 /i -1 Mpc) and 
with a sign predicted by theoretical models. This correlation could cause the existing 
lensing surveys at z ~ 1 to underestimate the linear amplitude of fluctuations by as 
much as 20 per cent depending on the source sample used, while for surveys at z ~ 0.5 
the underestimation may reach 30 per cent. The GI contamination is dominated by the 
brightest galaxies, possibly due to the alignment of brightest cluster galaxies (BCGs) 
with the cluster ellipticity due to anisotropic infall along filaments, although other 
sources of contamination cannot be excluded at this point. We propose that cosmic 
shear surveys should consider rejection of BCGs from their source catalogs as a test for 
GI contamination. Future high precision weak lensing surveys must develop methods 
to search for and remove this contamination if they are to achieve their promise. 

Key words: cosmology: observations - gravitational lensing - large-scale structure 
of Universe. 
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Weak gravitational lensing of distant galaxies has 
emerged as a powerful method for directly mea- 
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suring 



Melliei 
2003; 



the matter distribution in the uni v erse (e.g. 
119991: iBartelmann fc Schneider! 1200 it iRefregieJ 

following detection of the two-point func- 



tion of the lensing-induced shear by several groups 


iBacon et alj|200fj 


fVan Waerbeke et al. 2000: 


Rhodes et alJ 


20011: iHoekstra et alJ 12001 


Van Waerbeke et al.l 120021: 


Brown et al.l 120031 IJarvis et al 


2003). Also, weak lensing 



has been s hown to be a very promising tool for precision 
cosmology (Hu 2002; Hutcrcr 2002; Abazaiian & Dodelson 



2003; Bcnabcd & Van Waerbeke 2004: Bernstein & Jain 



2004 



Tereno et al] 12003: lUpadhve et al] 12003) . and has the 



clear advantage of being the simplest of the low-redshift 
cosmological probes to understand theoretically, since it is 
sensitive primarily to the dark matter whose interactions 
are described purely by gravity. Thus, most of the potential 
systematic errors in weak lensing are observational, i.e. 
associated with attempts to measure a small signal, in 
contrast with galaxy surveys, where large amounts of 
information are "lost" on quasilinear or nonlinear scales 
because of the lack of a robust prediction. 

Examples of the observational systematics in weak lens- 
ing include uncertainties in the point-spread function (PSF) 
of the telescope, errors in correcting the ellipticity of a 
source galaxy for the PSF effects, star-galaxy separation, 
selection biases, deblending errors, detector nonlinearities, 
and noise-related biases. The list is long, and much of 
the effort by the weak lensing observers in the past sev- 
eral years has been devoted to these issu es jKaiseill2000t 



Erben et all200lt [ Bernstein fc JarTis!l2002l:lHirata fc Seliakl 
2003al: IVale et al]l2004t IVan Waerbeke et al.ll2005l) . Nev~ 

theless, it is important to remember that there are astro- 
physical uncertainties associated with weak lensing. Some 
of these, such as the error in TV-body simulations or bary- 
onic cooling effects, are limited to small scales only and 
are likely to be redu ced significantly as the simulations 
improve dWhitdEiol . Another astrophysical uncertainty 
is the redshift distribution dN/dz of the source galaxies. 
At least at the brighter magnitud es, these can be deter- 
mined via spectroscopic redshif ts (iBernstein fc Jainlliool 
llshak fc Hiratal 120051: iMa et al.l 120051) . and at all magni- 
tudes the distributions can be tested by comparing the 
alaxy-galaxy lensing signal among different so urce samples 
Mandelb aum et al.l2005ai iHuterer et al l2005h . Future low- 
frequency radio surveys may also open up the possibility of 
obtainin g redshifts on many galaxies from the Hi 21 cm 
line fe.g. iRawlings et al]|2004T) . In short, while the redshift 
distribution problem is not solved, there is no fundamental 
impediment to an accurate determination of dN/dz. 

The final major astrophysical contaminant of weak lens- 
ing may be intrinsic alignments, i.e. correlations of the 
galaxy ellipticites with each other or with the density field. 
These correlations violate the usual assumption in lensing 
shear surveys that the source galaxies are randomly oriented 
so that any correlation between the apparent ellipticities of 
distinct objects is due to lensing. The purpose of this paper 
is to provide an observational constraint on these types of 
intrinsic alignments, and then to assess the implied contam- 
ination of upcoming cosmic shear surveys. 

There are two types of intrinsic alignments that can 
contaminate the cosmic shear power spectrum, namely the 
intrinsic ellipticity-intrinsic ellipticity (II) correlation and 



the gravitational shear-intrinsic ellipticity (GI) correlation. 
The II contamination is the easier to understand: it results 
from the possibility that two source galaxies are physically 
near each other and have correlated intrinsic ellipticities. 
This effect contributes an additive bias to any two-point 
shear statis tic, though as pointed out b v lKing fc Schneider] 
(2002) andEei mans fc Heavens! ]2003h . its effects can be 
minimized by using photometric redshifts to cross-correlate 
the shapes of galaxies at different redshifts. The GI contam- 
ination results from the possibility that, given two source 
galaxies, one is in front of the other. In this case, the in- 
trinsic ellipticity of the nearby galaxy may be correlated 
with the density field that lenses the more distant galaxy, 
thus yiel ding a spurious contribu tion to the shear two-point 
function llHirata fc S cliak 200l])- This contribution cannot 
be eliminated by selecting galaxies at different redshifts as 
for the II contamination, and in fact doing so will tend to 
increase any GI contamination. 

So far, the main methods used to estimate the in- 
trinsic alignments and assess their implications for cosmic 
shear have been theoretical results (analytical or simulation- 
based) and observations of low-redshift galaxies for which 
cosmic shear is negligible and any observed correlations must 
be intrinsic. Both methods have their limitations: the theory 
of galaxy alignments is subject to all of the uncertainties in- 
volved in the theory of galaxy formation, whereas the obser- 
vations must be extrapolated from observable redshifts z ~ 
0.1 to the redshifts z ~ 1 of the source galaxies used for lens- 
ing. Nevertheless, the theoretical pr edictions vary by more 
than an order of magnitude for II iCroft fc Metzlerll2000l: 
[ Heavens et al.l l2000l: ftee fc. Per] l2000l. |200ll : ICatelan et al] 
l200ll:ICrittenden et al.ll200ll : |jindl2002Tl . so there is a clear 
role for observations in distinguishing which of these models 
is correct. The theory of GI correlations is ev en more rudi- 
mentary: only simple analy tical models exist llHui fc Zhand 
l200aiHirata fc Seliaklliooil . 

At present, the most statistically powerful con- 
straint on II correlations co mes from the SuperCOS- 
MOS data ferown et al.ll2002i) . which were re-analyzed by 
iHevmans et al] ]2004f) with the conclusion that the II corre- 
lations are smaller than most of the theoretical predictions. 
The Sloan Digital Sky Survey (SDSS) is ideally suited to im- 
proving upon SuperCOSMOS because of the availability of 
spectroscopic redshifts, which allows for measurement of the 
correlations only between galaxies that are near each other 
in three-dimensional space, without the "noise" introduced 
by pairs of physically unassociated galaxies that happen to 
lie along the same line of sight. One can also attempt to 
constrain the GI correlation using low-redshift surveys such 
as SDSS. The GI correlation is essentially a measure of the 
correlation between the intrinsic ellipticities and the matter 
density; on sufficiently large scales, one can use the correla- 
tion between the intrinsic ellipticities and the galaxy density 
as a proxy. In this case the measurement of the GI corre- 
lation is model-dependent; however, we will argue that on 
large scales the associated uncertainty is probably less than 
the uncertainty associated with extrapolation to high red- 
shift. The ability of spectroscopic redshifts to isolate pairs 
of galaxies at the same redshift is valuable for GI correla- 
tion studies to reduce noise, just as for II. It also enables us 
to cleanly separate GI from the "spurious" (for this work!) 
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density-ellipticity correlations caused by galaxy-galaxy lens- 
ing. 

We note that there have been several observational 
studies of intrinsic alignments that were not motivated 
by potential contamination of cosmic she ar m easurements. 
These include bernstein fc Norberd ll2002l) and lHirata et all 
(2004), who were interested in intrin sic alignme n t con - 
tamination of g alaxy- galaxy lensing: [ Pen et alJ i2000^ . 
iLee fc Pen! lj200lfl . and ILee fc Penl (120021) . who were inter- 
ested in using the tidal influence on intrinsic alignments 
to rec onstruct the matter density field; and iNavarro et alJ 
(2004), who were interested in the formation of disk galax- 
ies. These results are difficult (or, in some cases, impos- 
sible) to interpret in the context of cosmic shear because 
they are made at too small a physical scale, use measures 
of the galaxy shape not easily related to ellipticity, mea- 
sure higher-order moments instead of density-ellipticity or 
ellipticity-ellipticity two-point functions, have complicated 
selection criteria that differ dramatically from those relevant 
to cosmic shear, or do not provide sufficient characterization 
of statistical errors. 

The outline of the paper is as follows. We begin by in- 
troducing the formalism used in this paper in ^21 The SDSS 
data used for this analysis is described in JO We describe the 
calculation of the relevant correlation functions and show 
the resulting signal in §4\ In we use the computed sig- 
nal to derive estimates of the contamination due to intrinsic 
alignments in current and future surveys. After proposing an 
explanation for the detected signal in SjU we conclude with a 
discussion of the implications of our results, and suggestions 
for future work. 

Here we note the cosmological model and units adopted 
for this work. Pair separations are measured in comoving 
h~ Mpc, with the angular diameter distance computed in a 
flat ACDM cosmology with Q. m = 0.3, Qa = 0.7. 



2 FORMALISM 

The form alism for the analysis o f the lensing shear two-point 
function dMiralda-Escude1ll99lD and the intrinsic alignment 
contamination is well developed. We will briefly summarize 
the important equations here, and then define some new 
variables that relate to observables in g alaxy surv e ys. O ur 
notation is consistent with that of lHirata fc Seliakl (l2004h . 

The observed shear 7 of a galaxy is a sum of two com- 
ponents: the gravitational lensing-induced shear 7°, and the 
"intrinsic shear" 7 7 , which includes any non- lensing shear, 
typically due to local tidal fields. Therefore the -E-mode 
shear power spectrum between any two redshift bins a and 
f3 is the sum of the gravitational lensing power spectrum 
(GG), the intrinsic-intrinsic, and the gravitational-intrinsic 
terms, 



CrM) = Cf E ' GG (a/3) + C7r- iJ M)+^ E '^(a/3). (I) 

The pure gravitational lensing term is given by the Limber 
integral 



iEE,GI/ 



C L {a/3) = / 



w -<x>w,bd ^ ^ (2) 

sm K x 

where \ is the comoving distance of the lensing screen, P& is 



the matter power spectrum, and the modified trigonometric 
functions are defined by 



flosinh^j, n K >0,R = 1 ^= 

sin K x = \ X, tt K = 

Ro sin , Q K <0,R = ? 



(3) 



and analogously for cosk %■ The window function for red- 
shift bin a is 



W a ( X ) 



-f2 m i?o(l + z) sin|f x 



x J fc (x) (cot if x ~ cot K x) dx'- 



(4) 



Here f a (x') = dn a /dx' is the distribution of comoving 
source distances for source sample a, and is equivalent to 
the redshift distribution if one fixes the homogeneous cos- 
mology x( z )- The intrinsic alignment power spectrum is 



P?i E (L csc K X, X) dx- (5) 



Here we have introduced the density- weighted intrinsic shear 
7 7 = (1 + 8 g )'j I , where S g is the galaxy (not matter!) over- 
density Pgj~p g — 1, and its i?-mode power spectrum P?i E ■ 
The GI or "interference" term is 

w a (x)Mx) 



c 



EE,GI 



(af3) = J 



P s ^i(Lcsckx) dx 



(6) 



To first order in the deflection angle, there is no B-mode 
induced by gravitational lensing shear. In this case, the B- 
mode power spectrum contains only a contribution from the 
intrinsic alignments, 

Cf s (a/3) 



fc,U)f(3(x) 

sinix 



P?i B {Lcse K x, x)dx- (7) 



There are lensing-induced contributions to the B-mode 
shear from other e ffects, such as multiple deflections 
llGoorav fc H J 120021: iHirata fc Seliakl l2003bi ICoorav et all 
2005) and, on small scales, modulat ion of the effectiv e source 
reds hift by galaxy cl u sterin g dSchneider et al.|l2002l) . 

IHirata fc Seliakl (l2004h provided formulas relating the 
power spectra in Eqs. l(2J _ Q to the real-space correlation 
functions. These formulas will be extremely useful because 
we directly measure the correlation functions in SDSS. 1 If 
one chooses any two points in the SDSS survey, their separa- 
tion in redshift space can then be identified by the transverse 
separation r p and the radial redshift space separation n. 2 
The + and x components of the shear are measured with 
respect to the axis connecting the two galaxies (i.e. positive 
+ shear is radial, whereas negative + shear is tangential). 
Then one can write the density-intrinsic shear correlation as 



1 We chose to measure the correlation function rather than the 
power spectrum simply because this involved minimal modifica- 
tion of pre-existing and well-tested code. 

2 The redshift space separation is frequently denoted tt; we use II 
to avoid confusion since the number it = 3.14... appears frequently 
in this paper. 
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Pg&(k) = -2tt J ^ s+ (r p ,U)J 2 (kr p ) r p dr p dU, (8) 

where £s+(r p ,Tl) is the correlation function between S and 
74.. The intrinsic-intrinsic correlations are 

(*) = J[Z++(r p ,n)J+(kr p )+t xx (r p ,n)J-(kr p )} 

(9) 



p EE 



and 

P*i B {k) = j[Z ++ (r p ,Tl)J_(kr p )+Z xx (r p ,Tl)J + (kr p )] 

x 27rr p dr p dn, (10) 

where J±(x) = [Jo(x)±J4(x)]/2, and £++ and£ X x represent 
correlation functions of <yj. The existence of the fl integral 
suggests the definitions 

/+00 
&+(r„II)dII, (11) 
-00 

and similarly for w++ and w xx . 
3 DATA 

The data used he re are obtained from the SDSS. The SDSS 
llYork et alJ I2006T) is an ongoing survey to image approx- 
imately tv steradians of the sky, and follow up approxi- 
mately one million of the detected obj ects spectroscopically 



matcly one million ol the detected objects spectroscopically 
jEisenstein et al]l200ll IStrauss et alTl2002t iRichards etall 
l2002Th The imaging is carried out by drift-sc a nning 
the sky in pho tometric conditions jHogg et al.1 l200lt 
[lvezi^et_alj |2004|), in fi ve bands (ugriz) llFukugita et alJ 
Il996l: ISmith et alJ 12002ft using a specially designed wide- 
field camera Tckinn^t alJ Il998t) . These imaging data are 
the source of the LSS sample that we use in this pa- 
per. In addition, objects are targeted for spectroscopy us- 
ing these data iBlanton et alJ l2003al) 



and 



are observed 
same telescope 



with a 640 - fiber spectrograph on the 
llGunn et alJ 120051) . All of these data are processed by 
completely automated pipelines that detect and mea- 
sure photometric properti es of objects, and astrometri- 
calibrate the data (iLupton et alJ l200lt iPier et ail 
iTucker et aT] l2005ll . The SDSS is well underway,_ 
has had five major data releases jStoughton et al 





Abazaiian et al 
Adclman-McC arthv et alJ 



Finkbcincr et al 
This analysis uses 
the spectroscopically observed galaxies in th e Value-Added 
Galax y Catalog, LSS sample 14 (VAGC; IBlanton et alJ 
2005), comprising 3423 square degrees with SDSS spectro- 
scopic coverage. 



same as t he L3-L6 su bsamples used in the weak lensing 
analysis of lSeliak et alJ <l2005al) . and are very similar to the 
L4-L 7 subsamples us e d in t he galaxy power spectrum analy- 
sis of iTegmark et all (l2004t) . 3 The absolute magnitude cuts 
are defined in terms of h = Ho/ (100 km s _ Mpc -1 ) such 
that one can implement the cuts without specifying the 
value of Ho- The magnitudes have been corrected for extinc- 
tion, K- correct ion, and e voluti on. The extinction correction 
uses the ISchleeeT et al.l <ll998T) dust redde ning maps, with 
the ex tinction-to-reddening ratios given in IStoughton et ail 
(2002). The if-correction to z = 0.1 is perfor med using the 
KCORRE CT vl_ll software as described bv IBlanton et alJ 
(2003b). The luminosities are corrected for passive evolu- 
ti on to z = 0.1 assum ing constant dM r /d,z = —1.6, as 
in IBlanton et al ] l|2003alh We exclude galaxies lying inside 
the bright star mask. Random catalogs were generated tak- 
ing into account the variation of spectroscopic complete- 
ness with position. The random points were assigned abso- 
lute magnitudes drawn from a distribution consistent with 
the real sample, and random redshifts were assigned for a 
given magnitude give n the selection function <j>(M\z) from 
IBlanton et al ] J2003d) . The sample selection criteria in this 
work and in that one are sufficiently similar that when 
4>(M\z) from that work is used to generate random redshifts, 
their probability distribution matches the redshift distribu- 
tion of the real lenses used in this work (within the noise). 
Shape measurements were obtained via re-Gaussianization 
for 96 per cent of this sample (see H3.2H : the remainder failed 
due various problems, such as interpolated or saturated cen- 
ters. 



3.2 Ellipticity data 

In addition to a sample of galaxies, we also need their el- 
lipticities. For this purpos e, we use the measurements by 
iMandelbaum et al.l (l2005alh who obtained shapes for more 



than 30 million galaxies in the SDSS imaging data down to 
magnitude r = 21.8 (i.e. far fainter than the spectroscopic 
limit of the SDSS). This secti on briefly describes the RE - 
GLENS pipeline presented in IMandelbaum et al.l (l2005af) . 
and the one modification we made for this paper. 

The REGLENS pipeline obtains galax y images in the r 
and i filters from the SDSS "atlas images" (IStoughton et al] 
2002). The basic principle of shear measurement using these 
images is to fit a Gaussian profile with elliptical isophotes 
to the image, and define the components of the ellipticity 

_ 1 - (b/a) 2 



(e+,< 



1 + (b/a) 



: (cos 2(f>, sin 2<J>), 



(12) 



where b/a is the axis ratio and <j> is the position angle of 
the major axis. The ellipticity is then an estimator for the 
shear, 



3.1 Galaxy subsamples 

In order to determine the correlation functions £s+, £++, 
and £ x x , one needs a sample of galaxies with which to mea- 
sure the intrinsic shear, and a sample of galaxies with which 
to trace the density field. For this paper, we use only the 
SDSS spectroscopic galaxies, divided further into luminos- 
ity subsamples by absolute magnitude in the r filter. The 
four subsamples are described in Table These are the 



(7+»7x) = 27£<( e +! e x)}> 



(13) 



3 The lTegmark et alJ 120041) sample differs in that a redshift cut 
was imposed to make t he sample volume-limi ted, which simplified 
their analysis. Also, the lTegmark et al.U2004l) subsample number- 
ing is different: their L4 is similar to our L 3, their L5 to our L4, 
etc. Fin ally, both the lTegmark et aJ] (2004) and the Sclia k et"all 
12005a) samples included a slightly stricter apparent magnitude 
cut than the sample in this work. 



© 0000 RAS, MNRAS 000, 000-000 



Intrinsic alignments in SDSS 5 



Table 1. The luminosity su bsamplcs used in th i s analy sis. The number of galaxies in each sample does not include failures. The linear bias 
is from lTegmark et alJ 120041) and lSeliak et all j2Q2£a|), normalized to erg = 0.9. The "total" linear bias and (2) are the number-weighted 
averages of the other samples. 



Subsample 


Magnitude 




Number of 


Number of 


Linear 


(*) 






range 




galaxies 


failures 


bias 




L3 


-20 ^ M r 


+ 5 log 10 h < 


-19 


66 312 


2 723 


0.85 


0.07 


L4 


-21 < M r 


+ 5 log 10 h < 


-20 


118 618 


4 601 


0.94 


0.11 


L5 


-22 < M r 


+ 5 log 10 h < 


-21 


73 041 


2 829 


1.08 


0.15 


L6 


-23 < M r 


+ 5 log 10 h < 


-22 


7 937 


307 


1.59 


0.21 


Total 








265 908 


10 460 


0.98 


0.12 



where 1Z ~ 0.87 is called the "shear responsivity" and repre- 
sents the resp onse of the ellipticity (Eq . 11211 to a small shear 
ijKaiser et alJll99.4 iBernstein fc ,Tarvisll2002fl . In practice, a 
number of corrections need to be applied to obtain the el- 
lipticity. The most important of these is the correction for 
the sme aring and circularization of the galactic images by 
the PSF; lMandelbaum et all <2005al) uses the PSF maps ob- 
taine d from stellar images by the psp pipeline llLupton et alJ 
l200lf) . and c orrects for these using t he re-Gaussianization 
technique of iHirata fe Seliakl i2003al) . which includes cor- 
rections for non-Gaussianity of both the galaxy profile and 
the PSF. A smaller correction is for the optical distortions 
in the telescope: ideally the mapping from the sky to the 
CCD is shape- preserving (conformal), but in reality this is 
not the case, resulting in a nonzero "camera shear." In the 
SDSS, this is a small effect (of order 0.1 per cent) which 
can be identified and removed using the astrometric solu- 
tion JPier et all2003ft . Finally, a variety of systematics tests 
are necessary to determine that the shear responsivity 1Z 
has in fact be en determined correctly. W e refer the inter- 
ested reader tofMandclb aum et alJ i2005aT) for the details of 
these corrections and tests. 



On e modification to the IMandelbaum et alJ 
d2005aT) pipeline is neces sary for the analysis here. 
Mandclbaum ct al rejected all galaxies with 

photo's CR or INTERP flags set, i.e. if photo had to 
interpolate over a cosmic ray, bad CCD column, or bleed 
trail in order to obtain the galaxy image. These cuts 
are reasonable for the faint galaxies used as sources in 
galaxy-galaxy lensing; however, the spectroscopic galaxies 
typically have large area and high probability of having 
one of these defects present somewhere within the image. 
Indeed, 33 p er cent of the spectrosco pic galaxies are 
rejected by the IMandelbaum et al] i2005al) cuts, and of this 
subset, 88 per cent come from CR and INTERP. Moreover, 
signal-to-noise is not an issue for the spectroscopic galaxies, 
with typical detections of ~ 100er (versus ~ 10-20cr for the 
photometric galaxies at the limit r — 21.8), so one does 
not worry about noise-related problems when doing the 
interpolation. We therefore replaced this cut with a cut on 
INTERP .CENTER, eliminating only those for which the 
defect was very close to the object's center and therefore 
may have had a relatively larger effect on the determination 
of the centroid or other parameters. 



4 COMPUTATION OF CORRELATION 
FUNCTIONS 

4.1 Estimator 

Our estimator for the the intrinsic alignment correlation 
function is a generalization of the LS (|Landv fc Szalavll993h 
estimator for the galaxy correlation function. The LS esti- 
mator for the galaxy correlation function £(r p , II) is 



(D - Rf _DD~ 2DR + RR 



RR 



RR 



(14) 



where DD is the number of pairs of real ("data") galaxies 
with separations r v and LT, RR is the number of such pairs 
in a random catalog, and DR is the number of pairs of data 
and random galaxies with this separation. (DR and RR are 
understood to be rescaled in proportion to JV g ai and iVg al 
respectively if the number of random catalog galaxies differs 
from the number of data galaxies, which is usually the case.) 
From our perspective, the key advantage of the random cat- 
alog subtraction in Eq. 11411 is that (D — R) = 0. This means 
that if there is any additive systematic error 8D in the data, 
we find that the bias in the correlation function is 



<<S£(r p ,II)> = 



2(D - R)5D + 5D 2 
RR 



(15) 



since (D — R) = 0, the bias is thus second-order in any sys- 
tematic. This feature makes the LS estimator more robust 
than previous estimators such as DD/ RR — 1, and is re- 
tained in most of the recent estimators for the correlation 
function or power spectrum. 

When we compute the galaxy-intrinsic shear correlation 
function £ g +(r p , TV), we can generalize Eq. 11411 to 

I ( m S + {D-R) S+D - S+R 
£ 9+ (r Pi n) = — = — , (16) 

where S+D is the sum over all pairs with separations r p and 
II of the + component of shear: 



S+D : 



E 



2TZ 



(17) 



where e+(j\i) is the + component of the ellipticity of galaxy 
j measured relative to the direction to galaxy i, and 1Z is 
the shear responsivity. S+R is defined by a similar equation. 
Averaged over a statistical ensemble, (S+) = (D — R) — 0, 
so the cancellation of systematics to first order (Eq. 1151 
also applies here. The subtraction of S+R in Eq. llbl is 
identical with the usual random catalo g subtraction proce- 
dure in galaxy-galaxy lensing studies (ISheldon et al 



proce- 
j 12004 
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iMandelbaum et al"ll2005alb We emphasize that positive £ g+ 
indicates a tendency to point towards overdensities of galax- 
ies (i.e., radial alignment, the opposite of the convention in 
galaxy-galaxy lensing that positive shear indicates tangen- 
tial alignment). 

For the intrinsic shear-shear correlation functions 
£++(r p , II) and £ x x (r p , IT), we simply use the estimators 



Galaxy-intrinsic shear correlation function 



where 



S+S- 



S+S+ 
RR 



and £xx = 



RR 



= E 



e+(j\i)e+(i\j) 
(2ft) 2 



(18) 



(19) 



and similarly for S X S X - Since (S+) = (S x ) = 0, the can- 
cellation of systematics to first order works again, i.e. the 
square of any spurious source of shear adds to Eq. 118H in- 
stead of the shear itself. 



4.2 Implementation 

We wrote two pipelines to compute the correlation func- 
tions, described below. Pipeline I is based on a tree correla- 
tion function code kindly provided by Ryan Scranton, and 
expanded to include quantities with spin (e.g. ellipticities), 
the use of bins in II in addition to r p , and to make the 
jackknife subregions have sides that are roughly equal; the 
jackknife method is used to obtain error estimates. Pipeline 
II is a "brute force" correlation function code designed to be 
as simple as possible; its main purpose is to provide a com- 
pletely independent check on the much more sophisticated 
Pipeline I. Pipeline I was used for the main scientific results 
of this paper (except where indicated), but we have checked 
that these are not significantly altered by using Pipeline II 
instead (see the result section). 

In order to find pairs of galaxies, Pipeline I uses the 
SDSSpix package. 4 To avoid noise in the determination of 
D — R in the S+ (D — R)/ RR estimator of the GI correlation 
function, we use 10 random points for each real galaxy in 
the catalog. The correlation functions are computed over a 
120 /i _1 Mpc (comoving) range along the line of sight from 
n = -60 /i _1 Mpc to n = +60 /i -1 Mpc, divided into 30 bins 
4 /i _1 Mpc in size, and the projected correlation function is 
computed by "integration" (technically summation of the 
correlation function multiplied by AIT) over IT. This calcu- 
lation is done in Nun = 10 radial bins from 0.3 < r p < 60 
h~ 1 Mpc. (Note, however, that to avoid calculating pairs over 
excessively large separations, we imposed a cut such that the 
maximum angular separation considered is equal to that at 
60 /i Mpc at z = 0.05, so the relatively small number of 
lower-redshift pairs do not contribute to the result at the 
largest values of r p .) Covariance matrices are determined 
using a jackknife with 50 regions. This number was cho- 
sen to be large enough to obtain a stable covariance matrix 
for the fits (it must be larger than JViV 2 ; see Appendix D 
of iHirata e t al. 2004}) but small enough that the size of a 
given jackknife region is larger than the scale on which the 



URL: http : //lahmu.phyast .pitt . edu/~scranton/SDSSPix/ 
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Figure 1. The correlation functions w g +(r p ), w++(r p ), and 
nixx(rp) obtained from the L3, L4, L5, L6, and the full galaxy 
samples. Each of the 10 bins contains the same range in r p for 
each of the samples, but some of the error bars have been slightly 
displaced horizontally for clarity (L5 has not been displaced). The 
L6 data have been multiplied by 0.1 so that they can fit on the 
same scale. All the errors are 68 per cent confidence bands, and 
the errorbars are highly correlated on large scales. 



correlation is to be measured. Our results for the galaxy- 
galaxy correlation function ^(r p ,Il) and the projected cor- 
relation func t ion as a function of luminosity match those in 
IZehavi et alJ ((2005), a similar analysis of SDSS data cover- 
ing a smaller area of the sky. 

Pipeline II is simpler in that for density-shape correla- 
tions, the S+D estimator is used instead of S+D~ S+R; and 
that there is only one bin in II, All = 120 ft - 1 Mpc. Because 
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Figure 2. The correlation functions w g +(r p ), u>+-|_(r p ), and 
«i x x (fp) obtained from the L3— L6 galaxy samples. Same as Fig. HI 
except that Pipeline II was used, and results are not shown for 
the full sample. 



this pipeline does not use D — R in any of its estimators, its 
calculations were completed with only one random point for 
every real galaxy instead of 10. 



4.3 Results 

The correlation functions obtained via Pipeline I are shown 
in Fig. and via Pipeline II in Fig. [5] We note that the dif- 
ference between the two results is significantly less than lcr, 
indicating that these two independently written pipelines 
give the same result. In order to c onvert w ox into A^ /j n t 
in g-g lensing for comparison against Irlirata et alJ J2004h , it 



is necessary to multiply by -1 (to account for the different 
ellipticity convention) and divide by All + w gg (r p ), where 
AIT = 120/i -1 Mpc an d the projected corr elation function 
can be estimated from IZehavi et al.l ^2p05) , for which the 
luminosity subsamples have n early the sam e selec tion as in 
this work. See Appendix A of iHirata et alJ ll2004l) for more 
details. Note that the results are still not quite ap propriate 
to compare against those in IHirata et al.l ll2004l) , because 
the Aji n t so obtained from, e.g., our L6 results here would 
be that expected if one did g-g lensing with L6 lenses and 
sources only, instead of L6 lenses and all sources. 



4.4 Systematics tests 

4-4-1 45- degree tests 

The 45-degree rotation of all sources is one of the standard 
tests of systematics. We compute both w gx and w+ x in a 
manner analogous to the computation of the signal described 
in H4.2I using Pipeline I only. This computation was done for 
all four luminosity bins. We defer analysis of the results of 
these systematics tests to H5.1. II where we show consistency 
of the 45-degree signals with zero. 



4-4-2 Large line-of-sight separation 

Another test that was performed was to compute the signal 
in the usual manner, but instead of integrating along the 
line of sight from —60 to 60 ft _1 Mpc, we integrate in the 
two ranges satisfying 30 ^ |II| < 90 h~ Mpc. If any sig- 
nal is found for w g +, w++, or w X x when computed in the 
usual manner, computing it over large line-of-sight separa- 
tions with a null result can allow us to verify that the signal 
is truly of astrophysical origin rather than due to some sys- 
tematic. Results described below (" H5.1.1H reveal no evidence 
of a systematic. 



4-4-3 Jackknife errors 

In this section, we consider the accuracy of the jackknife 
errorbars, a crucial part of determining the significance of 
any detections and of placing constraints on non-detections. 
There are two competing concerns in determining the num- 
ber of jackknife regions: the first is that there should be suffi- 
cient regions , where Nun = 10) that the covariance 
matrix is not too noisy; the second is that the regions must 
be large enough that they are truly statistically independent. 
If the regions are too small compared to typical pair separa- 
tions, then the lack of statistical independence could lead to 
underestimation of errorbars, and therefore overestimation 
of the significance of detections or overly tight constraints 
on non-detected quantities. 

In order to test whether our errors are overly influenced 
by edge effects, we redid the comparisons of w g +(r p ) (us- 
ing the S+D/RR estimator), w++(r p ), and w X x{r p ) for L3 
(which is most prone to edge effects due to its low redshift) 
using 25, 50, 75, 100, and 200 jackknife regions. Our expec- 
tation is that if edge effects are truly significant, then we 
should see the size of the errors at large r p decrease when 
the number of jackknife regions is increased; our hope is that 
the errors will have converged for our choice of 50 regions. 
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Error on w (r ) relative to 50 jackknife regions 
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Figure 3. The size of errors of w g + for L3 as a function of r v for 
the density-shape correlations, relative to the size of the errors 
for 50 jackknife regions. 



We also note that jackknife covariance matrices indi- 
cate a high correlation between radial bins on large scales, 
where the correlations are higher for the lower luminosity 
(and therefore lower redshift) bins since a given transverse 
separation corresponds to a larger angular separation. For 
example, the correlation coefficient for w g + between the two 
outermost radial bins is 0.90, 0.95, 0.73, and 0.45 for L3, L4, 
L5, and L6 respectively. All analysis in the following sections 
includes the full covariance matrices so that the correlations 
are taken into account. 

As shown in Fig. [3] while the errors in w g +(r p ) are 
clearly noisy at the 15 per cent level (which is unsurprising, 
since 50 is not that much larger than 10 3 ^ 2 w 32), there do 
not appear to be any significant edge effects. While the er- 
rors for the largest numbers of samples do clearly decrease 
on large scales relative to the errors with smaller numbers 
of samples, it is not apparent that there is any statistical 
significance to this decrease, and even if it is significant, 
it is not a problem for our choice of 50 jackknife regions. 
We have also confirmed that edge effects are not important 
for w++(r p ) and Wxx(r P ). We can thus conclude that our 
choice of 50 jackknife regions should not cause any system- 
atic underestimate of the errors. We address the question of 
whether our use of 50 jackknife regions may lead to overly 
noisy covariance matrices that affect our results in .5.1.11 



5 ESTIMATES OF CONTAMINATION 

In this section we estimate limits on the intrinsic alignment 
contamination for the current and future cosmic shear sur- 
veys. Our computations here use the vanilla ACDM cosmol- 
ogy with the parameters obtained bv lSeliak et alJ <l2005bl) : 
fib = 0.046, fi m = 0.28, n s = 0.98, as = 0.9, and 
Ho — 71.0 km s _1 Mpc -1 . Note however that the fractional 
error bars on our intrinsic alignment measurements are much 



greater than those on the cosmological parameters, so rea- 
sonable variations in the cosmology should not substantially 
affect our results. We begin by fitting several intrinsic align- 
ment models to the data ( H5.ll . We then consider methods of 
extrapolating these models from the SDSS redshift z ~ 0.1 
to the redshifts z ~ 1 of interest for cosmic shear ( H5.2H . 
Then the implications for current ( .5.31 and future ( .5.41 
cosmic shear surveys are discussed. 

5.1 Model fits 

We consider several models for the intrinsic alignments, 
which can be used to estimate the amount of contamina- 
tion to weak lensing power spectra. The following types of 
models are considered: 

1. Power law fits: The simplest model is a power-law fit to 
the intrinsic alignment correlation functions w g +, w++, and 
Wxx- 

2. HRH* model: This model is designed to make o ur re- 
sults comparable to those of iHevmans et all ^2004), who 
pr esented this modified v ersion of the HRH model originated 
bv lHeavens et alJ teOOOT) . It fits Eg. I2T1 below to the shape- 
shape correlation measurements with a free amplitude. 

In any of these models, once w g +, and w X x are 

determined, we can convert to P g ~i, P? r E , and P?j B using 
the Hankel transform relations (Eqs. In the latter two 

cases, the power spectra are immediately useful for intrinsic 
alignment studies. In the former case, a conversion from g 
to 5 must be applied. The simplest method here is the linear 
bias assumption, 8 = g/b g , which is valid f or some bias b a 
on sufficiently large scale s. The values of b g llTeemark et alJ 
12004, ISeliak et alll2005a_ are given in Table and we at- 
tempt to determine the scales on which t he linear bias as- 
sumpt ion is trustworthy using figure 9 in iTasitsiomi et all 
(2004), which shows the stochasticity and bias as a func- 
tion of scale determined in terms of correlation functions, 
i.e. r = £ sm / '^/.__£mm and b 2 = £ 99 /£ mm , determined us- 
ing N-body simulations. As shown there, the stochasticity 
is consistent with 1 for all scales larger than 400 h~ 1 kpc. 
The bottom panel indicates that the bias is approximately 
linear for scales larger than about 5 ft -1 Mpc, but even at 1 
h~ Mpc scale, the bias is only about 25 per cent different 
from its value in the linear regime. 



5.1.1 Power law 

The power law approach is to fit the signal to the equation: 



w a+( r p) = A 9+ 



lh- 1 Mpc 



(20) 



and similarly for u>++ and w x x . The fits were done by com- 
puting the x 2 on a grid in (A g + 1 a g +) using the full jack- 
knife covariance matrix; because this covariance matrix is 
noisy, the usual \ 2 analysis must be modified. Usually, it 
is assumed that A\ 2 (i.e. the difference in x 2 between the 
best-fit values of the parameters and the true values) has a 
X 2 distribution with a number of degrees of freedom equal 
to the number of parameters being fit. This is no longer 
th e case with a nois y covariance matrix. In Appendix D 
of iHirata et alJ i2004h . we examined the distribution of the 
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Table 2. Best-fit parameters for power-law fits A[r p /(lMpc/h)] a 
to the intrinsic alignment signal; the 95 per cent confidence inter- 
vals shown here may include no constraint on a if the amplitude 
is consistent with zero at this level. 



Sample 


function 


A (h~ 1 Mpc) 


a 


L3 


w g+ (r p ) 
w ++ (r p ) 

W V V (Vr>) 


0.0821°;™ 

-0.018+°:°^ 
0.005+°,° 3 ° 

0.022 


-0.181- 
-1.131- 
-0.68+°° 

— OO 


LI 


w g+ (r p ) 
w ++ (r p ) 
Wxxt'p) 


0.020±g;JSg 
(-5.71H2) x 10-5 
(3.8ll 5 |) x 10- 4 


0.0131- 
-5.5l- 

_7 1+5.8 


L5 


w g+ (r p ) 
w++(r P ) 
i»xx (r P ) 


o.ao±g;g 

n +0.035 
U.UJ1_q Q31 

n ml +0.030 

u.un_ 029 


-0 66+ ' 54 
-1 9+ 12 
-2.4+°° 

— OO 


L6 


w g +(r P ) 
w++(r p ) 
i»xx(r P ) 


Q Q+3.5 

o.o 2 2 
04+° ;45 
-0 25+ 1 ' 05 


n 77 +0.29 
u - ' '-0.30 

-l-8l~ 


All 


w g+ (r p ) 
w ++ (r p ) 
»xx (r P ) 


n nq«+ - 067 
(4.3 + | j) x 10" 3 
(7.2+ 7 3 2°) x lO" 3 


n 5Q+0.65 
_ 2g + oo 
O 1 +oo 

z * A -oo 



jackknife or bootstrap-derived x 2 f° r correlation functions 
computed in JV bins and M regions; a similar argument al- 
lows one to determine Ax 2 for TV bins, M regions, and K 
parameters. In our case, TV = 10, M = 50, and K = 2, and 
the 75th, 95th, and 99th percentiles of Ax 2 correspond to 
Ax 2 = 4.20, 9.55, and 15.44 instead of 2.77, 5.99, and 9.21 as 
derived from the standard x 2 distribution. Thus, our confi- 
dence regions are larger than those that might be predicted 
naively using the standard \ 2 distribution. The resulting 
fits and Ax 2 contours (which are only strictly accurate for 
Gaussian error distributions) are shown in Fig. 01 and the 
best-fit parameters are given in Table [5] with 95 per cent 
confidence intervals derived using Ax 2 = 6.17 (the value 
for one variable constraints). We defer further discussion of 
the apparent detection of nonzero w g +(r p ) in L5 and L6 to 
section |S| The last row in Fig. [I] shows confidence regions 
derived from the full sample of galaxies. It is clear that the 
majority of our constraints come from L4. 

We also performed the power-law analysis on the 45- 
degree rotated signals, w gx (r p ) and w+ x (r p ), for each lu- 
minosity bin. Of the 8 signals, 7 were consistent with zero 
at the 68 per cent confidence level for —4 ^ a ^ 1, and 1 
was inconsistent with zero at that level but was consistent 
with zero at the 95 per cent confidence level. Consequently, 
we conclude that this test does not reveal any systematics 
contaminating our analysis. 

Finally, we performed the power-law analysis for 
w g +(r p ) in L5 and L6 (the two bins with apparent detec- 
tions) with large line-of-sight pair separations only, as de- 
scribed in H4.4.2I The signal computed in this manner is 
consistent with zero at the 68 per cent confidence level for 
— 4 sC a ^ 1, which confirms that the detected signal when 
computed in the usual manner is indeed of astrophysical 
origin. 

To demonstrate explicitly the goodness-of-fit of the 
power-law model, and the inconsistency with zero for L5, 



L6, and the full sample, Table [3] includes the best-fit pa- 
rameters, their \ 2 and p(< x 2 )-value, and Ax 2 = X 2 (^ = 
0) - x 2 (best fit) and its p(< Ax 2 ). (Note that p(% 2 ) and 
p(Ax 2 ) represent the probability for a random vector to 
yield lower values of x 2 or Ax 2 by chance, and as such is 
technically equal to 1 — p, where p is the probability that 
the x 2 w iU exceed the given value. Thus the p(< Ax 2 ) for 
the A = case versus for the null hypothesis tells us the 
confidence level at which our results exclude the null hy- 
pothesis.) Because of the noisiness of jackknife covariance 
matrices, the x 2 values do not actually follow the x 2 distri- 
bution, and the p(< \ 2 ) and p(< A% 2 ) values are therefore 
computed using a simulation that acco unts for this fact (fo r 
more information, see Appendix D of iHirata et al.ll2004t) . 
Also, to demonstrate that our results are not overly influ- 
enced by noisy covariance matrices due to the use of only 
50 jackknife regions, we have repeated the analysis for 100 
jackknife regions, and included the same information for that 
analysis in Table [3] Results for L6 use the distributions for 
9 radial bins; all other bins use 10. We note that since the 
use of 100 jackknife regions involves different weighting of 
the data, the actual best-fit parameters, Ax 2 values, and p- 
values are not expected to be precisely the same as with 50 
regions; however, if the fit is good, the best-fit values should 
be consistent with each other and the p- values for a robust 
detection should be small in either case. 

There are a few conclusions that we can draw from Ta- 
ble [3 First, any changes to actual best-fit parameters due 
to noisiness of the jackknife covariance matrix (which would 
tend to affect the results with 50 jackknife regions more than 
with 100 jackknife regions) are much smaller than the vari- 
ation allowed within the 95 per cent confidence intervals on 
those parameters. In both cases (50 and 100 jackknife re- 
gions), the x 2 values for all samples except L6 are on the 
small side, as reflected in the p(Ax 2 ) values tending to all 
be relatively small. This may reflect some overestimation of 
the errors and therefore underestimation of detection signif- 
icance. The larger x 2 for L6 may indicate some deviation 
from power-law behavior, but the actual value of x 2 and 
p(< x 2 ) is not sufficiently large that using this model is un- 
reasonable. Finally, we consider the Ax 2 values for the spe- 
cial case of A = relative to the best-fit x 2 - As shown, the 
results for L3 and L4 are consistent with zero, with the Ax 2 
not being significant. For both L5 and L6, high-significance 
detections were obtained whether we consider the Ax 2 val- 
ues using 50 or 100 jackknife regions. For the full sample, the 
detection significance is roughly 99 per cent in both cases 
(slightly higher for 50 regions, slightly lower for 100 regions). 



5.1.2 HRH* model 

The HRH* model llHevmans et all2004l) treats the intrinsic- 
intrinsic term according to: 



w ++ (r p ) 



where r = 



wxx(r p ) 
A 



1+1 - 

ro 



1 + {r/B) 



■dru, (21) 



+ r p . iHevmans et alJ l|2004h fit this model 

with the parameters 7 = 1.8, B — 1 h^ 1 Mpc, and ro = 
5.25 h~ Mpc, leaving A as a free parameter. Their results 
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Power law fits to intrinsic alignment signal 
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Figure 4. Power law fits to the intrinsic alignment signal from Pipeline I. From top to bottom, the rows represent the luminosity ranges 
L3, L4, L5, L6, and the full sample; from left to right, the columns represent the w g +, and w x x correlations respectively. The dots 

indicate the \ 2 minimum, and the contours represent 75, 95, and 99 per cent confidence regions. In the L6 density-shape plot, there was 
an insufficient number of pairs in the innermost bin to establish a reliable jackknife error estimate, so this plot is based on only N = 9 
data points. 



are shown in Table 0] one can see that a detection was 
achieved using SuperCOSMOS. The original Super COS- 
MOS intrinsic alignments analysis llBrown et al.ll20o3l used 
a median redshift of z me d = 0.1 for their sample with limit- 
ing magnitude b j = 20.5, and this numb er was also used in 
the re-analysis bv lHevmans et al.l II2004J). The source of the 
quoted z me( j is a parametric model from lr3augh fc Efstathioul 
( 1993) using APM data, which was obtained using fits to a 
particular form of the redshift distribution including a pa- 
rameter z me d(bj). It is possible that the weak dependence 
of Zmed on bj for bright galaxies assumed by the fitting for- 



mula may have unduly influenced the fits, causing 2 me£ j to 
have been underestimated at intermediate magnitudes. 

However, more recent data from the 2 degree Field 
(2dF) galaxy survey has enabled a direct measurement 
of the bj band lumi nosity function and k + e-corrections 
llNorberg et alj [2002). from which we find z me d ~ 0.16 for 
this limiting magnitude. In support of this calculation, we 
note that from the 2dF catalog itself (which has a lim- 
iting magnitude of bj = 19.45), z me d for bj < 18.5 is 
0.082, and for bj < 19.45 is 0.11. If we extrapolate the 
bj < 19.45 number to the SuperCOSMOS flux limit using 
the Euclidean method, z med ~ l0°- 2b - r i "™ (valid for z <C 1) 
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Table 3. Best-fit parameters for power-law fits A[r p /(lMpc//i)] a to the GI intrinsic alignment signal u> 9 + (r p ) (with 95 per cent confidence 
intervals), the x 2 an d P(< X 2 )> an d Ax 2 = X 2 (A = 0) — x 2 (best fit) and its p(< Ax 2 ). These numbers are all shown for both 50 and 
100 jackknife regions. 



Sample 50 jackknife regions 100 jackknife regions 







best-fit 
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0)-(best fit) 




best-fit 
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0)-(best fit) 




A (h- l Mpc) 


a 


x 2 


p(< X 2 ) 


AX 2 


p(< Ax 2 ) 


A (ft-iJVIpc) 


a 


x 2 


p(x 2 ) 


Ax 2 


p(< Ax 2 ) 


L3 


0.082i°;^ 9 


-0.18±~ 


3.9 


0.09 


9.6 


0.95 


062 +0 ' 157 
u.uoz_ Q 12 i 


-0.08l^ 


1.2 


0.14 


3.9 


0.80 


LI 


020"* -0 ' 115 
u.uzu_ og5 


0.013+^ 


5.7 


0.23 


2.1 


0.51 
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u.uzo_ Q 102 
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, we obtain z me d = 0.174. If we determine an exponential 
Zmed ~ 10 a bj > lim and extrapolate that to bj t i im = 20.5, we 
obtain z me d = 0.144. Both of these extrapolations merely 
serve as sanity checks of our result from the luminosity func- 
tion calculation, and the fact that z met j for bjj im = 19.45 
is already larger than 0.10 in the 2dF data further supports 
our assertion that the SuperCOSMOS calculations signifi- 
cantly underestimated 2 me d- The large underestimate of the 
median redshift for the SuperCOSMOS data increases the 
theoretical predictions for the observed signal, and hence 
decreases the b e st-fit value of A; we therefore conclude that 
iHevmans et alJ ll2004h underestimated both the value of A 
and its error bar. We have re-computed the theoretical pre- 
dictions for the ellipticity variance at z me d — 0.16 and found 
the correction to be a factor of 3.0-3.4 over the range of an- 
gular scales considered by SuperCOSMOS (25-93 arc min- 
utes). The value shown in the table has therefore had a re- 
scaling factor of 3.2 applied (using 3.0 or 3.4 instead changes 
the results by < 0.2cr). Note that the change in the red- 
shift distribution does not affect the statistical significance 
of the SuperCOSMOS detection, although there are poten- 
tial sources of systematic error in this data set that could 
lead to a spurious signal. 

The constraint on A from the full SDSS sample is 
slightly tighter than that from SuperCOSMOS, but due to 
its lower central value does not provide a detection. Combin- 
ing the two results is dangerous because the southern stripes 
of SDSS overlap with the SuperCOSMOS survey area and 
hence the two results are not truly independent. The results 
in Table El suggest that the SDSS and SuperCOSMOS con- 
straints on A are mutually consistent, with any discrepancy 
below the la level and thus not very statistically significant, 
particularly in light of the large systematic uncertainty due 
to the uncertainty in median redshift. Indeed, this compari- 
son highlights the need for accurate redshift distributions in 
order to properly constrain the amplitude of intrinsic align- 
ment correlation functions. 



5.2 Redshift evolution 

In order to estimate the amount of intrinsic alignment con- 
tamination in cosmic shear surveys, one must extrapolate 
from the SDSS redshift range z ~ 0.05-0.2 to the redshifts 
of the sources used in these surveys, z ~ 0.3-2. We attempt 
several types of extrapolation here. Note that in each case 
the extrapolation becomes successively more dangerous as 
we go to higher redshift: our results for the lowest-redshift 



Table 4. The intrinsic alignment amplitud e A in Eq . 1211 . as de- 
termined in this paper and bv lHevmans et all 120041) . All results 
are 95 per cent confidence; the error bars from this work are de- 
termined using the jackknife confidence intervals as described in 
J8-l.il The SuperCOSMOS and CO MBO-17 ellipticity data are 
described in lBacon et alJ J200fj) and lBro wii^^dj j2003j ), r espec- 
tively. The analysis of SuperCOSMOS bv lHevmans et all i2004l) 
gave a value of A = (0.9 ± 0.5) X 10 -3 ; we have re-scaled this 
result by a factor of 3.2 to account for the more recent determi- 
nation of the SuperCOSMOS redshift distribution. See i|5.1.2l for 
details. 



Sample 


A/10" 3 


Reference 




SDSS L3 


-1.0 ±6.0 


This work 




SDSS L4 


1.7±4.3 


This work 




SDSS L5 


5.8 ±7.4 


This work 




SDSS L6 


-43 ± 122 


This work 




SDSS L3-L6 


1.8 ±2.3 


This work 




COMBO-17 


< 5.4 


Hevmans et 


al. (2004) 


SuperCOSMOS 


2.9 ± 1.6 


Hcvnians et 


al. (2004) 



surveys such as CTIO are probably very good, and any at- 
tempt to make statements at redshifts z ~ 2 are probably 
close to meaningless. 

Extrapolation in redshift involves three major issues. 
One is that the intrinsic ellipticities of galaxies may change 
with time due to (for example) mergers. A second issue is 
that galaxies move, and thus the intrinsic alignment power 
spectra can change even if the orientation of each galaxy 
remains fixed. Finally, we have established that the intrinsic 
alignment signal depends on the sample of galaxies consid- 
ered, and it is not clear which sample of galaxies observed at 
2 = 0.1 is most similar to (or is evolved from) a particular 
sample of galaxies at higher redshift. Each of these issues 
must be addressed in the context of both the density-shape 
and shape-shape correlations. 

In t he case of the shape-s hape correlation, the argu- 
ment in iHirata fc Seliakj (120041) would suggest that in the 
linear regime, w ++ (r p ) and w X x(r P ) would be redshift- 
independent (again assuming the intrinsic ellipticities of in- 
dividual galaxies do not change). Of course, most of the 
constraints from lensing, and most of the observational data 
presented here, are from nonlinear scales where the galaxy 
correlation function £ (r p ,II) ^ 1. In the nonlinear regime, 
IHevmans et alJ i2004ll argued that the correlation function 
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of the unweighted intrinsic shear f 1 would vary slowly with 
redshift, while the density-weighted correlation functions 
w++(r p ) and iii xx (rp) would increase at low 2 due to the 
growth of perturbations 5 g . In particular, if one examines 
the strongly nonlinear regime £ ss ^> 1, and assumes stable 
clustering with £ 99 oc (1 + z)~' i+1 (7 ~ 1.8 is the slope of 
the galaxy correlation function), then it follows that both 
w++(r p ) and w X x(r p ) should scale as (1 + z)~' i+1 . 

However, if one takes seriously the concept of stable 
clustering, then the unweighted intrinsic shear correlation 
function, , should be constant when measured in phys- 

ical rather than comoving separation. In this case, the cor- 
relation functions of 'y 1 at redshift z are 

£ ++ (rv,rn; Z ) = (l + ^r 3 £ ++ f^,JI_;o) . (22) 
Integration over ry then gives 

w ++ (r p ;z) = (l + z)- 2 w ++ (yTI ;0 ) ' (23) 

and similarly for w x x (?>)■ 

In the case of the density- shape correlation , the sim plest 
prescription was provided by lllirata fc Seliakl (|2004) , who 
argued that if the intrinsic ellipticities of individual galaxies 
did not change, then on linear scales ws+(r p ) would scale as 
the growth factor. On the other hand, on nonlinear scales, 
stable clustering would suggest 

ws+(r P ;z) = (l + z)- 2 w s+ (yT7 ;0 ) ' (24) 

in analogy to Eq. il'i.'it . The analyses below consider both 
the linear evolution and the stable clustering methods for 
extrapolating the signal to high redshift. We do not consider 
other methods of extrapolation since these would likely be 
much more complicated and it is not clear that they would 
have higher fidelity. 

5.3 Current surveys 

In the past several years, a number of cosmic shear sur- 
veys have reported results for erg, which is the cosmological 
parameter most readily accessible to cosmic shear. It is ap- 
parent that there is some disa greement between th e deter- 
minations (see, e.g., Table 1 of lHevmans et al1l2005l) . While 
the differences may arise due to shear calibration biases of 
up to 15 p er cent between the PS F co rrection methods used 
(see, e.g.. lHirata fc Seliakll2"003al and lHevmans et alJl2005l) 
or uncertainties in z m , there is also the possibility that the 
samples used for each of these studies were contaminated by 
intrinsic alignments in very different ways. Unfortunately, 
the determination of how our results apply to these sam- 
ples is highly nontrivial. We have selected our galaxies in 
the r band at z ~ 0.1, but the cosmic shear selection at 
2 ~ 0.6-1.0 is usually in the observed-frame red/far-red 
bands (similar to SDSS r or i), which is frequently in the 
rest-frame ultraviolet. This makes it very difficult to de- 
termine whether the samples in those tables correspond to 
the lower-contamination L < L* samples, or the higher- 
contamination L > L* samples. 

In principle, one possible method of determining the 
correspondence between samples would be to use a simple 



one-parameter family of templat es from the stellar popula- 
tion synthesis code described in iBruzual fc Charlod (l2003ll 
with exponential star-formation rates. Even this simple 
model can be shown to accurately reproduce the colors of 
SDSS main spectroscopic sample galaxies over for < 2 < 
0.25. One could then find the regions in the color-color and 
color-magnitude diagrams corresponding to each of the lu- 
minosity bins, use the models to evolve them in redshift, 
and see what they looked like at 2 ~ 0.6-1.0. However, such 
a procedure entails enough systematic uncertainty that it 
is unclear whether the results would be trustworthy. Hence, 
we defer a detailed discussion of how the detected GI align- 
ments contaminate the results of current surveys to future 
work. 

It is interesting to note that on small scales the GI 
effect on the observed angular shear correlation function 
may have been detected in the R ed-sequence Cluster Sur- 
vey fRCS). iHoekstra et alJ IJ2004) observed an anticorrela- 
tion between the shapes of bright (19 < Rc < 21) and 
faint (21.5 < Rc < 24) galaxies. This anticorrelation is not 
consistent with either cos mic shear or estimates of the sys- 
tematic errors in RCS, and lHoekstra et alJ i2004f) attributed 
it to the alignment of the bright galaxies with the elliptic- 
ities of their dark matter haloes (which lens the more dis- 
tant, faint galaxies). 5 The halo ellipticity will tend to lead 
to a stronger lensing signal along the lens major axis, i.e. 
stronger source tangential ellipticity along this axis, and 
therefore an anticorrelation of lens and source ellipticities. 
This halo ellipticity effect is a special case of the density- 
intrinsic shape correlation, although the physical explana- 
tion in terms of alignment of galaxy light with halo mass 
distribution suggests it would be limited to scales below the 
virial radius of the halos. On larger scales, such an anti- 
correlation can be interpreted as ellipticity of local large- 
scale structure rather than of individual ha los. Of course, 
the galaxy samples in IHoekstra et alJ i2004l) were selected 
by putting most of the "lenses" at lowe r redshift than the 
"sources," which Ulirata fc Seliakl (2004) argued maximizes 
the GI effect. Nevertheless, precisely this type of separation 
of source screens has been proposed for cosmic shear tomog- 
raphy studies. The apparent detection of halo ellipticity also 
suggests that in some cases, the GI effect can exceed 100 per 
cent of the GG effect and result in negative shear correla- 
tion functions. Of course, such regimes must be avoided for 
cosmic shear studies. 



5.4 Future surveys 

The GI contamination estimated using the power-law fits is 
shown in Fig. |S] and the II contamination in Fig. HJ for a 
cosmic shear survey with a median redshift of z me d = 1.0. 
The full sample (L3-L6) GI results should be treated with 
caution because they are obtained by averaging the vari- 
ous galaxy samples whose GI contamination amplitudes are 



5 iMandelbaum et al did not find such a correlation in 

the SDS S, but the samples of galax ies used were very different, 
and both IMandelbau m et all 1200 5b') and this work suggest that 
these correlations are strongly dependent on the sample of galax- 
ies considered. 
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not consistent with each other. The average is thus appro- 
priate to the SDSS spectroscopic sample, or other samples 
at the same redshift range with the same selection crite- 
ria, but may not be a good guide otherwise. In particular, 
when considering GI contamination for a future survey, one 
must determine which of the luminosity subsamples is most 
similar to the sources to be used in that survey. 

Fig. |7|shows the estimated fractional GI contamination 
for a survey with z met j — 0.6 (determined using linear evolu- 
tion, though both methods give nearly identical results for 
these low redshifts) only for those samples that had a sta- 
tistically significant determination of signal. A comparison 
between Figs.|^|and|7|shows that the contamination is more 
severe at lower redshifts, as expected. 

We note that, as discussed in H5.ll these estimates of 
contamination are not reliable on small scales due to the 
breakdown of the linear bias assumption. At < 1 /i~ 1 Mpc, 
or L >~ 1000 at z me d — 1, we may be overestimating the 
bias by approximately 30 per cent, and consequently under- 
estimating the contamination by that amount. 



5.5 Cosmic shear tomography 

Thus far we have only discussed the effects of intrinsic align- 
ments on cosmic shear autocorrelation studies. In this sec- 
tion, we also consider the effects on cosmic shear tomog- 
raphy, the cross-correlation of pairs widely separately in 
redshift, either via down- weighting pairs close in photomet- 
ric redshift, or by explicitly separating into multiple sam- 
ples based on photometric redshift or apparent magnitude 
such that the samples lie in distinct regions in redshift. In 
principle, while the cross-correlation should still give a GG 
signal due to lensing by structures between the observer 
and the lower-redshift source sample, any intrinsic align- 
ment II contami nation will be highly su ppressed. However, 
as pointed out in lHirata fc Seliakl (12004^ . the GI contamina- 
tion may still be significant when selecting pairs separated 
in redshift space, because the galaxy at lower redshift can 
be aligned with the same tidal field that lenses the higher- 
redshift galaxy of the pair. In light of our robust detection 
of the GI contamination using pairs at the same redshift, 
more investigation into the effect of this contamination on 
cosmic shear tomography is warranted; this investigation is 
beyond the scope of this paper. Some of these issues have 
been investigated bv lKinsf <l2005t) . 



6 INTRINSIC ALIGNMENTS AND BCGS 

An early study of alignments o f cluster g alaxy ellipticities 
with the cluster major axis, iBinggelil 1^982), found a 
tendency for brightest cluster galaxies (BCGs) to align 
with the cluster ellipticity using 44 Abell clusters with 
z < 0.1. L a ter st udies confirmed this "Binggeli effect," 
iFuller et alJ (1 19991) with poor clusters, IWest fc Blakeslee 
( 2000) in 3 dimensions for the Virgo cluster. and lKirri^t^n" 
(2002) using ~ 300 clusters in SDSS data over a large 
range in redshift, 0.04 < z < 0.5. This effect has been 
explai ned i n terms of anisotropic infall along filaments 
JWesdll989l) . an explanatio n that is support ed by numerous 
iV-body simulations (e.g. iDubinskv l Tl998l . [The related 
"Holmberg effect" iHolmbenir^oa) . the correlation of 
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Figure 5. The allowed range of GI contamination for each lu- 
minosity subsample. The power-law fits are used, with the left 
column showing the results for the stable clustering assumption 
and the ri ght column show i ng the results for linear evolution as 
argued by |Hirata fc SeljalJ i2004) . The bottom and top curves 
show the 95 per cent confidence region assuming a power law in- 
trinsic alignment model with index —3 < ct g + < +1. The center 
curve shows the contamination predicted by the best-fit parame- 
ters in Tablel^l The median source redshift assumed is 2 me( j = 1.0. 
[The constraints on ct 3 + are imposed because for « 9 + ^ +1 or 
a g + ^ —4, the Hankel transform defining P s ~i (k) becomes ill— 
defined. A cutoff value greater than —4 was chosen because oth- 
erwise the correlations at very small scales dominates the power 
spectrum. Note that in the cases of L5 and L6 where we have a 
detection, a g + is constrained to lie within the range given here 
at > 99 per cent confidence.] 
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satellite positions with the major axis of the primary 
(not in clusters), has also been studied f a irly extensively 
obser va tionally llHawley fc Peebles! Il975t IValtonen et alJ 
" iMacGillivrav et al.1 Il982t IZaritskv et al] Il997t 



Sales & Lamba 2004; Brainerd 2005) and with simulations 
jNavarro et al] 20041: lAzzaro et al.ll2005l: [Rang et aljl200ll 



iLibeskind et al J 120051: IZentner et al] 12003) . and is another 
possible explanation, though unlike for clusters the results 
- magnitude and even sign of the effect - have been con- 
flicting.] This effect may be used to explain the detection 
of nonzero w g +(r p ) in L5 and L6: if some of these bright 
galaxies are BCGs of clusters, then we expect to find an 
overdensity of galaxies along their major axis, leading to 
positive w g +. In support of this hypothesis, we note that 
the majority of these bright galaxies are red, and many pass 
the color cuts to be included in the SDSS spectroscopic 
LRG sample; it is well known that BCGs tend to be red 
galaxies. 

Naturally, BCG alignments can only explain a detection 
on small (<~ 2 h~ Mpc) scales. How ever, when combined 
with the result IIHopkins et al] 120051) that clusters them- 
selves are aligned with each other for separations up to f 00 
/i -1 Mpc, we can also explain the larger scale alignments, 
since it means that a BCG of a given cluster also tends 
to point preferentially towards overdensities of galaxies on 
larger scales and not just within its own cluster. We also 
note the result shown there (from iV-body simulations) that 
the alignment increases with increasing redshift, indicating 
that if our explanation of the effect is correct, then it may 
be worse for higher-redshift surveys than what has been de- 
tected here. 

If indeed the major source of density-shape correla- 
tion is BCG alignment, it may be advantageous for cos- 
mic shear surveys to reject BCGs, particularly in the lower- 
redshift bins where their intrinsic alignments can contam- 
inate the cross-correlation tomography signal. A variety of 
algori thms exist to identify clusters and their BCGs (e.g . 
iGladders fc Yeel200(ilGoto et al.l20ollBahcall et all2003) . 
and it is beyond the scope of this paper to examine how 
much the density-shape correlation can be suppressed by 
rejecting BCGs selected by each of these methods. 

IPereira fc Kuhnl ]2005l) use X-ray selected clusters with 
SDSS photometric and spectroscopic imaging data to find a 
tendency for radial alignment of cluster galaxies, not limited 
to BCGs. To explain their findings, they propose a paramet- 
ric resonance phenomenon that causes a tendency for many 
cluster galaxies to exhibit radial alignment relative to the 
cluster center. Their results suggest that removing BCGs 
from the sample may not be sufficient to eliminate GI con- 
tamination. Further investigation of our results, including 
a search for intrinsic alignments that distinguishes between 
satellites in clusters and the BCGs themselves, is necessary 
to understand the degree to which the effect they observe is 
important for cosmic shear surveys. 

We also note that the effect seen here is related to 
the effect seen by correlating positions of a fainter set 
of galaxies with the pos ition angles of lens galaxies in 
iMandelbaum et al] ]2005bl) . There, as shown in Figure 6 and 
the accompanying text, a statistically significant tendency 
for the fainter galaxies to be positioned along the major axis 
of the central galaxy was found for L5 and L6. 



7 CONCLUSIONS 

We have used the SDSS main spectroscopic sample to search 
for II and GI correlations by computing projected corre- 
lation functions w g+ (r p ), w ++ (r p ), and w xx (r p ) over the 
range of pair separations 0.3 < r p < 60 ft _1 Mpc. We have 
two main results coming out of these calculations. 

The first result is a constraint on the II correlation, 
expressed in terms of several models in H5.1 l and in terms of 
projected contamination for higher-redshift surveys in H5.4I 
As shown, we have no detection of II correlations either in 
luminosity bins or using the full sample. The limits rule out 
a significant contamination of shear signal by II, although 
contamination at the 10 per cent level for lensing surveys at 
z s ~ 1 is still a possibility. Thus, the II correlation may still 
be a concern for future cosmic shear surveys with expected 
statistical errors below the 1 per cent level, though it can be 
minimized by cross-correlating source galaxies in different 
redshift slices. 

The second result is a detection of GI correlations for 
L > L, galaxies (L5 and L6), as well as for the overall 
sample, with nonzero amplitude of the correlation function 
w g+( r p) f° r the power-law subscript —4 ^ a ^ 1, at the > 99 
per cent confidence level. We have no detection of GI cor- 
relations for L ^ L- t galaxies (L3 and L4), and have placed 
constraints on them for the first time. The agreement of our 
results from two independent pipelines to calculate the cor- 
relation functions confirms the robustness of our fundings. 
As for the II correlations, we have used several models of ex- 
trapolating these GI results to higher-redshift cosmic shear 
autocorrelation power spectra in i !5.4l to predict GI contam- 
ination at z s = 1 and L — 100 — 1000 of 5-25 per cent (L5), 
50-150 per cent (L6), or 0-15 per cent (full sample). These 
GI correlations are positive in sign, which is the opposite 
of the GG lensing induced signa l, so the interference is de- 
structive (Hirata & Scliak 2004). As a result the amplitude 
of fluctuations is systematically underestimated in current 
surveys where this effect has been ignored. The estimated 
error in the linear amplitude as from this effect is uncertain, 
but estimates presented here suggest current surveys under- 
estimate it by 0-20 per cent for z s ~ 1 and possibly up to 
30 per cent for shallower surveys with z s ~ 0.5. 

These results leave open a number of questions for fu- 
ture work. First, is there some way to minimize the contam- 
ination due to the GI correlations, which cannot be elim- 
inated by the scheme proposed for eliminating II correla- 
tions via cross-correlation of sources at different redshifts? 
We have proposed BCG alignments with cluster ellipticities 
as a possible explanation for the correlation, so eliminat- 
ing BCGs from the sample used for computing cosmic shear 
should be investigated as a means of reducing this align- 
ment, though as suggested bv lPereira fc Kuhnl (120051) . this 
step alone may not be sufficient. Second, what are the impli- 
cations of the GI detection for current cosmic shear studies, 
and could they explain the conflicting values of as (which do 
have other possible explanations, such as shear calibration 
bias or uncertainty in redshift distributions)? In order to 
investigate this question thoroughly, it will be necessary to 
determine the correspondence between the samples used for 
this work and for the current cosmic shear studies. Finally, 
it is clear that a systematic error at a 10 per cent level can- 
not be tolerated in future surveys such as LSST, PanStarrs 
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or JDEM, which require a sub-percent precision to achieve 
the stated goals. More detailed implications for future work, 
including cosmic shear tomography, which is susceptible to 
contamination from the GI detect ion, will be addressed in a 
following work. Ilshak et al.l (l2005Tl - 

In summary, our detection of the GI contamination of 
the cosmic shear power spectrum should serve as a useful 
starting point for further investigation into this effect, which 
is one of the main theoretical uncertainties in understanding 
current cosmic shear analyses, and may also lead to methods 
that can help reduce this contamination to the lowest level 
possible when interpreting data from future surveys. 
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Figure 6. The allowed range of II contamination for each lu- 
minosity subsample. The power-law fits are used, with the left 
column showing the results for the stable clustering assumption 
and the ri ght column show i ng the results for linear evolution as 
argued by [Hirata fc Seliakl i2004) . The bottom and top curves 
show the 95 per cent confidence region assuming a power law in- 
trinsic alignment model with index —1.5 < a < 0. The median 
source redshift assumed is z me d = 1-0- [The constraints on a are 
imposed because for «++ +0 or ct++ ^ —2, the Hankel trans- 
form defining P s -i (k) becomes ill-defined. A cutoff value greater 
than —2 was chosen because otherwise the correlations at very 
small scales dominates the power spectrum. We do not detect II 
in any of these cases, so the actual power law slope cannot be 
determined from the data.l 
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Figure 7. The allowed range of GI contamination for L5, L6, 
and the full sample determined using the power-law fits and lin- 
ear evolution for z me( ; = 0.6. The bottom and top curves show 
the 95 per cent confidence region assuming a power law intrinsic 
alignment model with index —3 < o 9 + < +1. The center curve 
shows the contamination predicted by the best-fit parameters in 
Table 
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